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, _ ' Abstract 

\o ■ 

^SJ ' Cosmic strings are considered in two types of gauged sigma models, which generalize the 

gravitating Abelian Higgs model. The two models differ by whether the U(l) kinetic term is of the 
Maxwell or Chern-Simons form. We obtain the self-duality conditions for a general two-dimensional 
^ I target space defined in terms of field dependent "dielectric functions". In particular, we analyze 

analytically and numerically the equations for the case of 0(3) models (two-sphere as target space), 
and find cosmic string solutions of several kinds as well as gravitating vortices. We classify the 
solutions by their flux and topological charge. We note an interesting connection between the 
^— J ■ Maxwell and Chern-Simons type models, which is responsible for simple relations between the 

' self-dual solutions of both types. There is however a significant difference between the two systems, 

, in that only the Chern-Simons type sigma model gives rise to spinning cosmic vortices. 

"o '. 

O^; PACS numbers: 11.27. +d, 04.20. Jh, ll.lO.Kk, ll.lO.Lm. 

^ '. 

; 1 Introduction 
• 1— I . 

' Topological defects are generally believed to have been formed during a series of phase transitions in 

^ , the early universe. In particular, many field-theoretical models suggest the formation of cosmic strings 

■ - - ' (for a review, see for instance Vilenkin and Shellard 1 ) . 

The BOOMERANG results '7, especially concerning the location of the acoustic peaks [HIE], essen- 
tially rule out GUT scale (~ lO^^GeV) topological defects being the sole source of the cosmic microwave 
background (CMB) anisotropy, which reflects the early universe density fluctuations, eventually leading 
to galaxy formation. However, recent analysis 0011113111 shows that a mixture of cosmic strings and 
inflation is consistent with current CMB data. 

Independently of a possible role in the large scale structure formation, and even if their mass scale is 
lower than the GUT scale, cosmic strings may have also a number of significant observable astrophysical 
effects. Just to mention a few: Cosmic strings could be sources of double images P, of gravitational 
waves ^1 ^3 and of ultra high energy cosmic rays . Recently it was even suggested ^3] that they 
can serve as gamma ray burst engines. Thus, cosmic strings are still of considerable interest in cosmology 
and astrophysics. 

The most popular field-theoretical system, used to model cosmic string generation 1[, is the Abelian 
Higgs model. There exist, however, other systems which may be used for the same purpose and are 
not more complicated. As long as we model cosmic strings as static cylindrically symmetric sources 
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coupled to gravity, which is essentially a two dimensional system, we may borrow circular solutions of 
other well-known two dimensional systems, namely non-linear sigma models. 

Non-linear sigma models |14| are very much used as effective theories describing various systems such 
as low energy effective QCD, and are also used as "toy models" for gauge theories. A generalization 
of these, where the global symmetry group or its subgroup is promoted to be local, are gauged sigma 
models which appear naturally in supersymmetric field theories. 

All these models have been extensively studied for decades. Most of the results were obtained in 
flat (Minkowski) spacetime, but more recently gravitating solutions were studied too. However, the 
gravitating sigma model solutions are usually taken to be spherically symmetric as in the cases of 
gravitating Skyrmions and textures (see e.g. Volkov and Gal'tsov JSl)) although some authors gave 
attention to cylindrically symmetric solutions of the non-gauged HTj as well as the gauged models 
with either Maxwell [HI or Chern-Simons terms [H EOl EH E2 ■ 

Here we will analyze the gravitating gauged sigma models with a general two-dimensional target 
space defined in terms of two field dependent " dielectric functions" , which may be viewed as generalized 
Abelian Higgs models: Maxwell type and Chern-Simons type, according to the gauge field term in the 
Lagrangian. We will find the general conditions for self-dual cosmic string solutions and get several 
kinds of interesting solutions. 

We will start by considering, in the next section, static and translationally invariant solutions of a 
generic gauged gravitating non-linear sigma model of a Maxwell type, and find simple conditions for the 
existence of self-dual solutions. In sections 3 and 4, we concentrate in the 0(3) model and add rotational 
symmetry to obtain self-dual cosmic string solutions. In section 5, we move to the Chern-Simons type 
of model, and consider stationary (not static) solutions. First we obtain the self-duality conditions for 
this kind of solutions, and then we specialize to spinning cosmic vortices. In section 6, we discuss the 
relation between the two types of theories and solutions, and we end with a discussion of the topological 
charges in section 7. 



2 The Generalized Abelian Higgs Model 

We start by considering a generalization of the Abelian Higgs model defined by the action: 

5 = 1 ds^ (^y^i\mD,'i>nD'^<f) - um) - y^rnDF^.F^-- + 



(2.1) 



where fid^j) and f2(|'I'|) are non- negative dimensionless functions, which may be interpreted as Weyl 
factors of a conformally-fiat target space metric of a non-linear sigma model with local U(l) symmetry 
|23j. The function i?2(|'i'|) plays further the role of a dielectric function j^EHI- In what follows, we 
use the collective notation £a, (a = 1,2), and refer to both as "dielectric functions". The action H2.1|l 
also generalizes the ones considered by Lohe [221 E|- Furthermore, the action 1)2. l|l arises in various 
extended supergravity theories 23: . In the general case, we may also add the term Ssd^De'^'''"^ Ff^i^Fp^j , 
but this term vanishes identically for the configurations considered in the present paper and has no 
effect on the field equations. 

The field equations derived from the action H2.1|) are: 



£U\^\)D D^^+ ^ D $n^$ + — + ^ F F^" = (22) 



{£2m)F''n = / = -'-e£,mm*iD-'^) - ^D''<i>r) (2.3) 

— i?^. + U^im) ((i?^$)*(Z?,<i>) + iD,^nD^<^>)) - C/(|$|)g^. + 



SttG 



£2m) { g'^^F^^F^, -f If^aF^V ) = (2.4) 



Conventions: Z?^ = - ieA^, signature (+, -, -) and R^^^ = d^Vl^ ~ d^J^I^ 
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Usually (but not always) we will be interested in potentials, which ensure spontaneous symmetry 
breaking and which lead to a massive gauge field. That is to say, potentials with a circle of degenerate 
minima at (say) |$| = u and with positive second derivative at v. We further normalize the potential 
such that it vanishes in the vacuum. Thus, altogether: 

U = 0, C/' = |||^0; for ^^v (2.5) 

With these normalizations, the Higgs mass and gauge boson mass are, respectively: 

, U"iv) 2 e^v'E^jv) 

In order to study cosmic string solutions, we assume the metric and matter fields to be static and 
symmetric under z-translations. Thus we assume that A^^ and $ depend only on the two transverse 
coordinates x'', and for the metric we use the following form: 

ds^ = N^{x'')dt^ - j^j{x'')dx'dx^ - K^{x^)dz'^ (2.7) 
We also require the presence of a magnetic field only, i.e. that the gauge potential will have the form: 

A^dx^" = A,{x^)dx^ (2.8) 

such that the Maxwell tensor will contain a single magnetic component B: 

F^^dxi" A dx" = ~By/yi\e,jdx'- A dx^ (2.9) 

where I7I = |det(7y)|. In order to get the field equations for static solutions, we compute also the 
components of the Ricci tensor: 

N 

Roo = —-\7,{KV'N) 
K 

R33 = ^V,{NV'K) 



R^0 = R^3=0, i?03=0 (2.10) 



Rr3 = i?y(7) + T7V.V,7V+-V,V,if 



where Vi is the covariant derivative with respect to the two-dimensional metric 7^ and Rij^j) the 
corresponding Ricci tensor. 

A significant simplification of this system is obtained if self-duality conditions are satisfied, i.e. if 
the system admits a Bogoniolnyi limit j2Hl ■ It is well-known that in the usual Higgs model, the flat 
space considerations can be carried over to curved background if iV(a;*) and if (a;*) are constants, say, 
1. We will see now that the present generalized Higgs model has also a Bogomolnyi limit if we keep 
N{x'^) = K{x^) — 1. If we use these conditions, we find that the (00) and (33) components of Einstein 
equations will be satisfied only if: 

um) = l£2m)B' (2.11) 

Now we turn to the (ij) components of the Einstein equations, or even better to Gy which vanish 
identically. Consequently, Tij = as well and we get: 

i£i(|$|)7'='(i^fc$)*(A$) + urn) - l£2m)B''^ 7., = (2.12) 

Using (|2.11|l . this condition simplifies further and it follows that it is equivalent to the curved spacetime 
version of the self-duality condition: 

- iv^/W\^^Jl"'Dk^ (2.13) 
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where rj = ±1 corresponds to self-dual or anti self-dual solutions. This is a first order equation for the 
Higgs field. Analogously, Ea. (|2.11|l is a first order equation for the gauge potential. 
We also find the following expression for the two dimensional Ricci scalar: 

i?(7) = -87rG(fi(|$|)7*^(A$)*(i?j$)+ 4(7(1$!)) (2.14) 

which serves as the Einstein equation for the two metric 7y . 

Equation 12.2|l for the Higgs field becomes a consistency condition, which constrains the form of the 
potential [/(|$|): 

,.|*|^,m)BH-i|i + f|| = ,2.15) 

Maxwell equations 1)2. 3|l give: 

d,iS2m)B) = -^£,m)dm' (2.i6) 

which is a second order equation for the gauge potential, but it is not an independent one, as is easily 
shown: If we think of i? as a function of |$|, we may use H2.11|l and get from (|2.15(l the following relation, 
which we will refer to as the " Bogomolnyi constraint" : 

(f2(|<i>|)i?)+77e|<i>|£i(|$|)=0 (2.17) 



This is actually again the Maxwell equation H2.16|) in disguise, and can be used together with (|2.11() to 
get the function J7(|<i>|) for any set of given "dielectric functions" £a, a = 1, 2. In order to do it directly, 
we need to express the magnetic field B in terms of the potential U. If B has a definite sign, we infer 
from the Bogomolnyi constraint that if2(|$|)-B(|$|) > for j] = +1, and similarly for negative values, 
so we may take a square root of (I2.11|) : 

B = T]^/2U/£2 (2.18) 

and get a simple relation between the potential and the two dielectric functions, which is a necessary 
condition for the action 1)2. l|l to have self-dual solutions, i.e. a Bogomolnyi limit: 



V2f2(|$|)J/(|$|)+e|$|£i(|$|) -0 (2.19) 



The simplest case where we can apply (|2.19|l is of course £a = 1, which reproduces immediately the 
standard Higgs potential: 

Um) = -y-\'^?f (2.20) 

where the Bogomolnyi relation between the coupling constants holds: 

a = eVa = 2 (2.21) 

Note that here, the vacuum expectation value v enters into the potential as an integration constant, 
and has also a role of the value of |<I>| for which the magnetic field vanishes. 

The next case is the non-linear 0(3) sigma model, which is obtained by taking Ei to be the usual 
conformal factor for with a radius /x/2 [Tlj : 

fi(|a>|) = l/(l + |<i>|V/i')' (2.22) 

The parameter sets a second energy scale in the system and renders £i dimensionless. It is very simple 
to integrate the Bogomolnyi constraint 1)2.17(1 . or equivalently Ea. H2.19|) . also in this case and get the 
following form of i3(|<I>|) and potential, which allows a Bogomolnyi limit in the 0(3) generalized Abelian 
Higgs model: 

""^i"^!^ - 2(1 + v^i„-^) s^m) 1 + 1*1 ^ ^ ^ 

Um)^ = ^- ' ' ' (2.24) 
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Here again v is an integration constant, which parametrizes the minimum of the potential and the field 
value for which the magnetic field vanishes. Generally there is no relation between this scale and the 
scale fi. We therefore have a one-parameter family of potentials in this 0(3) generalized Abelian Higgs 
model for any given f2(|*&|), which we may take at our will. Some special cases with £2(1*51) = 1 were 
already analyzed by several authors mainly in flat space [291 \'M)\ . but also with coupling to gravity |18| . 
A more "symmetric" picture may be obtained by using the angular variable on S'^ defined by: 

|<I>| = Aitan(e/2) (2.25) 

Generically, the potential exhibits a U(l)-symmetry breaking minimum at |<i>| = v, but there are two 
special limits v = 0, oo, where the ground state is only a point in target space (the north or south poles 
of S^); namely, no symmetry breaking occurs. This pattern of symmetry breaking is reflected by the 
Higgs and gauge boson masses, which turn out to have the following equal (due to the self-duality) 
values: 

cv 

rriH ^mA = — — tti-^ (2.26) 
f + I [l'^ 

The two cases of symmetric vacuum correspond to vanishing Higgs and gauge masses. 

Just as a final check we note that the Higgs potential f (|2.2U|l with (|2.21l) ') is obtained in the limit 
— > oo from (|2.24|) with £2 (1*1' I) = 1- The masses which are obtained in this limit are the usual ones. 

Now we return to the Einstein equation (|2.14|) and notice that for self-dual solutions, its right hand 
side is actually a two-dimensional divergence (it must be as we will see in section , and may be 
simplified to: 

^ levrG^^ (7^J£2(|$|)i?9, log |$|) (2.27) 

Without loss of generality, we may take the two-dimensional transverse part of the metric tensor to be 
conformally flat, i.e. 7ij = H^{x'')6ij so the Einstein equation H2.f 4|l reduces further to: 

S'^d, (^d,\ogH~^£2m)Bd,\og\'i>\^ = (2.28) 

This equation is not of first order, but it may get the form of a two dimensional Laplace equation, if we 
introduce a "super dielectric function" y^(|$|), which solves the following equation: 

dA 



l^l^=^^2(|*i>|)i?(|<i>|) = V2^(^M^ (2.29) 



or equivalcntly the second order equation: 



^ ^'*l:^)+^^^i(l*l)-0 (2.30) 



|$| V 

Now we can use the function A, in order to give Eq. (|2.28|) the form of a Laplace equation: 



[logH- — ^ ) = (2.31) 



Special solutions will be discussed in section 4. 



3 Flat Space Solutions 

First we discuss flat space solutions within this general framework. As mentioned above, some of the 
solutions are already known but new ones can be easily obtained. 

In order to study a single cosmic string solution, we take the usual cylindrically symmetric Nielsen- 
Olesen ansatz for n flux units: 

$ = ;i/(r)e"'^ , A^dxf ^ A{r)d'4} (3.f) 
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where is a second energy scale, which is generally independent of v. We will also assume that the 
dielectric functions depend on |<i>| through the dimensionless ratio |$|//i = / only. For further use we 
define here v/^i = (3. 

As for the boundary conditions, we will see that they will not be the same for all systems, but will 
rather have to be adapted to the specific system. However, since we are interested in solutions with finite 
energy per unit length and finite flux, the boundary conditions at infinity should ensure asymptotically 
vanishing energy density. The usual Nielsen-Olesen conditions should be generalized such that all three 
contributions (scalar gauged kinetic term, potential and Maxwell term) will vanish asymptotically. That 
is: 

lim t/(/(r)) = 0, lim £2ifir))B^{r) = 0, lim f{r)£i{f{r))ieA{r) ~ n)^ ^ (3.2) 

r — >QO 1 — >oo r — ^oo 

In order to proceed, we will concentrate in the 0(3) model. The field equations (with rj = +f which we 
will mostly use from now on) are the following simple first order set: 

^ = eA{r) - n (3.3) 



(3.4) 



2(l + /32)f2(/)(l + /2) 

where a prime denotes differentiation with respect to r. 

We still have a freedom in the function £2, but we choose until further notice £2{f) = 1- The simplest 
case, which yields the closest to the Abclian Higgs flux tube, is the case (3 — 1 which corresponds to a 
potential with a minimum along the equator of . We thus impose the following additional boundary 
conditions: 

/(0)-0, A(0) = (3.5) 
and the two last conditions at 1)3. 2|l may be replaced by: 

lim A{r) = n/e (3.6) 

r — ^CJO 

Note that the boundary conditions we took enforce n to be negative and A{r) to be a non-positive 
decreasing function. The magnetic field B{r) is non-negative and the flux (in units of 27r/e) is —eA{oo) = 
— n, which is positive. 

There is no known analytical solution of this system, but it is very easy to get numerical solutions. 
We will comment a little more about it at the end of this section. Fig^ contains the field variables (G, 
A and B) in this (3=1 case. This is actually the solution discussed already by Mukherjee [30]. There is 
another related solution, which is obtained by reflecting this solution in the target space with respect 
to the equator, i.e. starting at / = cxd on the z-axis (r — 0) and decreasing with r to / = 1: 

lim /(r) = 00, lim /(r) = 1 (3.7) 

r — >0 r — *cxj 

The reflection property is much more transparent in terms of the angular field variable O, Eg. 1)2.25(1 . in 
terms of which the boundary conditions on the scalar field are: 

9(0) = TT, lim e(r) = 7r/2 (3.8) 

r — ^00 

This reflected solution has obviously negative magnetic flux. 

Another solution, which was discussed already by Schroers corresponds to the limiting case 
of potential with a minimum which does not break the U(l) symmetry. This corresponds, in our 
terminology, to /3 = while /i stays finite so w = 0. In this case, the minimum of the potential is on 
the north pole of S^. Finite energy (per unit length) solutions still exist, but they have now different 
boundary conditions namely: 

6(0) =7r, yl(0)=0 (3.9) 
lim e(r) = , lim B{r) = (3.10) 

r — ^00 r — >oo 

The last boundary condition in 1)3. 2|) is automatically satisfied as / ^ asymptotically. This gives a 
solution, which maps the z-axis (r = 0) to the south pole of and the circle r ^ cxd to the north pole. 
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Since there is no symmetry breaking at spatial infinity, the magnetic flux is not quantized but rather 
takes continuous values. 

In order to compare these solutions with the /3 = 1 ones, it is more instructive to present the "anti- 
Schroers" solutions, which are solutions to the same problem with the reflected potential (/3 — > oo). The 
boundary conditions should be therefore also reflected: 

9(0) = 0, A{0) = (3.11) 

lim e(r) = TT , lim B{r) = (3.12) 

r — !-oo r — >oo 

where again the last boundary condition in (|3.2|) is automatically satisfled, as now f'^Si asymptot- 
ically. There are two unusual features of this family of solutions: First, as mentioned above, there is no 
flux quantization, and second, there does not exist a solution with |n| = 1 but with \n\ > 2 j29j . Con- 
sequently, the energy density is maximal not on the symmetry axis, but rather on a cylindrical surface 
whose radius is of the order of 1/e/i. Similar behavior is seen in the analogous |ri| > 2 Nielsen-Olesen 
flux tubes 1 . 

Another kind of interesting potential is a triple well potential, whose minima are at the north and 
south pole and at the equator of the target space. In term of the angular field 8, it has the following 
simple form: 

C/(e)-^sin2(2e) (3.13) 

This potential is obtained by using f2(|$|) = id^l/M + ^/l'^'l)^ and w = ^ in Ea. l2.24|l . Actually, £2 
may be multiplied by an arbitrary constant, thus having the height of the potential free. Since the 
"polar" minima (at Q — 0,tt) do not break the U(l) symmetry, we do not expect flux quantization for 
solutions that approach asymptotically one of these minima. On the other hand, solutions which have 
G = 7r/2 as their asymptotic value, will have of course quantized magnetic flux. 

There are thus several kinds of finite energy (per unit length) solutions in this system. For 0(0) = 0, 
all three kinds with asymptotic values 6 = 0, ■7r/2, tt are realized. The fields 9, A and B in figures 1^5171 
which actually correspond to the Chern-Simons type system analyzed in sec. [S] depict also the solutions 
of the present system with 9(oo) = 7r/2, tt, respectively. We will say more about it in sec. |B| 

Similar solutions exist with 8(0) = tt. A different kind of solution starts with a "non-polar" value 
(i.e. 8(0) 7^ 0, tt), but for continuity must have n = 0. The asymptotic values correspond to each of the 
three possible vacua. 

The numerical procedure we have used is essentially the same as in previous works |31l I32| , and 
we will not elaborate about it more than the following remark regarding the treatment of boundary 
conditions. A common practice in numerical analysis is to replace the boundary conditions at infinity 
with the same boundary conditions at some finite point, r^o- This works well on the system (|3.3f) - (l3.4f) . 
provided that /? is not or cxo. If we try to shoot from Too to in these two cases, it is easy to see (from 
the equations in terms of the angular field, 8), that the solutions become constant. We are therefore 
forced to replace the boundary condition on 8 with 8(cx)) = P ± 6 where (5 is a (positive) parameter, 
which is to be determined numerically. It turns out, that the continuum of solutions is found by giving 
S slightly different values. 

Another remark concerns the classification of solutions according to the possible choices of boundary 
conditions. As is already clear from the flat space solutions, the 0(3) system has a reflection symmetry 
with respect to the S^ equatorial plane (i.e. 8— >7r — 8or/— >1//) together with (3 1//?. We may 
therefore limit ourselves to solutions with /(O) = (which for rj ~ +1 are also n < solutions) and 
classify them according to the asymptotic values /(oo), which are the various minima of the potential 
function. The n > solutions for the same rj and /3 (and hence the same potential), which start at the 
"south pole" of target space, are easily obtained from the n < solutions with /3 — > 1//3 and further 
obvious changes as n ^ —n and A — A. 
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4 Gravitating Solutions 



In this section we consider gravitating 0(3) cosmic strings. It is in fact very easy to get the gravitating 
self-dual cylindrical solutions (i.e. single cosmic string solutions) in the Bogomolnyi limit, since the 
field equations for the matter (scalar and gauge) fields are very similar to those in flat spacetime. More 
precisely, they are given by the following first order equations: 

^ = r/ (eA(r) - n) (4.f) 

which is identical with the Minkowskian one (for 77 = +1 compare ()3.3|l 'l. and: 

A' 7/e^2 p2 _ J2 



mr 2(l + /?2)£2(/)(l + /2) 
Einstein equations for the metric field (|2.28(l get in the cylindrical case the following form: 



(4.2) 



(logi7)'-^i3£2(/)(log/)') 
■qe J 



-I I 



= (4.3) 



Since we are interested here in cases where the matter fields tend asymptotically to their vacuum 
values, the geometry of space will evidently be conic with a deficit angle given by the usual relation: 

^ = - lim(r(logi7)') (4.4) 

We may obtain more information about the deficit angle by integrating equation (|4.3|) once and using 
the result at both ends, r = and r 00: 

|^ = ^|nS(0)|£,(/(0)) (4.5) 

where we used the boundary conditions H3.5|l and 1)3. 6|l supplemented by H{0) = 1, H'{0) = 0, and the 
fact that —nB{0) is always (i.e. for 77 = ±1) positive. The factor i?(0)£2(/(0)) in the last equation can 
be expressed in terms of the dielectric function £i{f) using H2.17|l and assuming a vanishing magnetic 
field as / ^ /3: 

^£ ^ 87TG^l^\n\ 1^ f£i{f)df (4.6) 

For the 0(3) sigma model £i{f) under consideration here, this integral can be explicitly calculated, or 
we may rather use the already given expression of the magnetic field - Ea. (|4.2|l . Both ways we get: 

^ 4TTGv^\n\ 

Note that these two last equations are valid for the /(O) = solutions - see final note of previous section. 
The well-known usual Higgs result is obtained for (3 = 0, while holding v fixed. On the other hand, 
we may take another limit to obtain the deficit angle for the " anti-Schroers" solutions by /3 — > 00, while 
holding fi fixed. This gives {5lp/2tt)s = 47rG'/i^|n|. This result applies also to the Schroers solutions. 
In Fig|21 we show the fields (8, A, B, H) for the case /3 = 1, SnGji^ = 0.5; compare with FigQ] 



5 Chern-Simons Type of the Generalized Higgs Model 

The possibility of a second field dependent "dielectric function" £2, used at the end of section|21 is much 
more interesting from another aspect, which is a connection with D = Z Chern-Simons theory. It turns 
out that the self-dual solutions for a, D = i system of a gauged sigma model coupled to pure Chern- 
Simons theory, are related to those of our generalized Higgs system. In other words, we may replace the 
Maxwell term with a Chern-Simons term in the Z? = 3 version of the action H2.f |) to get a generalized 
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Chern-Simons-Higgs (GCSH) model. We will find interesting relations between the self-dual solutions of 
this GCSH system and the previous one, which from now on we call generalized Maxwell-Higgs (GMH) 
model. 

There is of course a physical difference between time-independent solutions in both models, which is 
the presence of electric charge and electric field in the GCSH model. However, due to the linear relation 
between the magnetic field and the electric charge density, a |<i>|-dependent magnetic energy density 
appears, which in flat space mimics exactly the Maxwellian energy density in the GMH system, as will 
be shown below. This sheds new light on the well-known self-dual solutions of the Chern-Simons gauged 
0(3) sigma model [^|211I2S1ISZI- The self-dual solutions of the Chern-Simons type of the Abelian Higgs 
system |SHlEn]j fit also to this framework by the obvious choice ^ri(|$|) = 1. In order to demonstrate 
these relations we write the D = 3 action: 

s^J d'xy^i (^^£,mKD,'!>nD^^) - um) + i^^^^a^'- + le^^) (^-^^ 

where all the geometrical terms here are 3-dimensional, and k is a positive parameter (it cannot be a 
field-dependent function, as in the GMH case, because it would break gauge invariance) with dimensions 
of 1/length. Actually, the fields here have different dimensionalities than in D = 4, but the self-dual 
solutions of the GMH system are identical in D = 3 and D = 4, so the comparison becomes trivial. 
The field equations of this system are: 

£,m)D,D''^ + 7^4^ ^M*^"* + 7^4^ = (5.2) 
2|<i>| |$| d|<i>| 

= r = -'M\m'^>*iD'''^>) - HD'^m (5.3) 
^i?^. + ^im) {{D^^nD,<^>) + (i^^$)) - 2um)g^, = (5.4) 

Spontaneous symmetry breaking occurs if the potential has a circle of degenerate minima, just as for 
the GMH system. The mass of the Higgs field is given by the same expression as for GMH - Ea. H2.6|l . 
but the mass (squared) of the gauge field is now ^: 

if^)' (5,5) 

In order to find static self-dual solutions in flat space, one may proceed by dropping the curvature 
terms and look for conditions for Ea. (|2.13|l to be satisfied. However, we prefer to take the more general 
approach of dealing with the gravitating system from the beginning. We therefore try to repeat what 
we did in the previous sections and require 500 = 1 in a D = 3 version of the static metric (|2.7|l . It 
turns out that this cannot be done " naively" , but some modifications are required. The reason is the 
well-known property of Chern-Simons vortices, which should carry also electric charge so the Maxwell 
tensor must contain in this case also an electric field, and the vortex carries angular momentum as 
well. This addition has important consequences if we treat the gravitational field as dynamical, since 
the angular momentum of the Chern-Simons field forces the gravitational field to be stationary and not 
simply static. We therefore parametrize the D — 3 metric by: 

ds^ = N^ix'') {dt + L,{x^)dx'f - -iij{x^)dx'dx' (5.6) 
and write for the gauge field: 

A^dx^' = Ao{x'')dt + A,{x'')dx' (5.7) 
In order to get the field equations for stationary solutions, we introduce also the following notation: 

A, = A, - AqL, 

Di = di~ ieAi 

Lij = diLj - djLi = \f\ri\e-ijl 

F,j = d,A,~d,A,^-^\e,,B (5.8) 



^Compare Deser and Yang | 
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and compute the components of the Ricci tensor: 

Roo = ~N\/,\/'N ~ —f 

= ^e'''3j(A^^^) (5.9) 

i?y = Iv^V^'TV + i(i?(7) - N'^e)-i'^ 

where Vi is the covariant derivative with respect to the two-dimensional metric ^ij and i?(7) the 
corresponding Ricci scalar. 

Now we impose N — 1 and find that the (00) component of Einstein equations (|5.4II will be satisfied 
only with the following condition on the potential: 

U{\^) = _ (5.10) 

From the {ij) Einstein equations (or better from the G'-' equations), we find the form of the self-duality 
equation in this case: 

A$ = i??/Me»n'''^fe* (5.11) 
as well as the following expression for the two dimensional Ricci scalar: 

i?(7) SttGs (fi(|$|)7'^(A$)*(A$) + 4t/(|$|)) (5.12) 

Due to the self-duality, the other field equations simplify as follows. The non-diagonal Einstein equations 
give: 

da = 167rG377e|$|fi(|$|)Ao9,|$| (5.13) 
The spatial Chern-Simons equations give: 

a,Ao = -^|$|£i(|cI>|)9,|$| (5.14) 

K 

which is used in the time component of the Chern-Simons equations to obtain an expression for B: 

B^ii+-m^£,m))Ao (5.15) 

K 

Finally, the equation for the Higgs field gives another functional relation between the potential and the 
other quantities: 

rITI fi 

rjeimimB + ^ - Y(^o)'^(|*|'f i) = (5.16) 
Actually it is not an independent equation, since by substitution of (|5.10(l in (|5.16() . one may get back 

The two differential equations for £{x'') and Aa{x'^) may be also converted into equations for ^d*!*!) 
and Ao(|$|): 

167rG3rye|$|£i(|$|)Ao, (5.17) 



dm 

^^-^I.|.r(|.|) (5.18) 

and we find a simple expression of £ in terms of Aq: 

i^8TTG3K{c- Al). (5.19) 

where c is an integration constant, which should be non-negative for solutions with finite angular mo- 
mentum (where £ vanishes asymptotically). 
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As in the gravitating GMH model, the right hand side of (|5.12(l is a two-dimensional divergence and 
we find: 

i?(7) = 16^G3«:V, ( —Y'Aod, log |$| + ] (5.20) 

If we use a conformally flat metric, H5.20|l can be rewritten as: 



S'^d, l^dj logiJ ;^^o^J 1*1 j = ^T^G3.KcHH (5.21) 

Now we consider some special cases that allow self-dual solutions. We will then concentrate in 
rotationally symmetric solutions, i.e. gravitating vortices. 

5.1 ^i($|) = l 

This is the Chern-Simons version of the usual Higgs model, which we may consider either in flat space 
or coupled to gravity. Integration of Eq. (|5.18|l is trivial and gives: 

^0 = ^{v^ - m (5.22) 

where v is an integration constant. In flat space t — Q and Ea. H5.10|) gives immediately the well-known 
|38[ I39| sixth order potential: 

Um) = ^^\'^?{v^-\'^?? (5.23) 
In the case of a gravitating system, we solve for ^(|$|) and find: 

tm)-'^-^{-~{v'-Wf) (5.24) 

ri 

where tr = 4K^c/e^. The potential is easily found to be: 

8k 



- ^ - 1*1')' - ^G, {a - [v' - \^^fY (5.25) 



which has an unbounded additional term, thus rendering the whole system possibly unstable. Notwith- 
standing the possible instability, this potential was discussed by several authors |19[ 121)1 1211 \TB and 
vortex solutions were also obtained. There is a simple and natural way to cure the ill-behaved potential, 
which is generalizing this system to the 0(3) sigma model. 

5.2 ^i(<l>|) = l/(l+|<l>|V/i2)2 

This is the Chern-Simons gauged 0(3) sigma model considered in flat background by some authors 
EHl ESI El ^D- Within the present framework the analysis is straightforward. First we integrate 
Eq.(E3HIl and find: 

In flat space £ — and Eq. H5.10|) gives immediately the following potential: 

f/(|$|)-^^^^i^— (5 27) 

Some authors have already discussed the self-dual solutions in flat space for special cases of the potential, 
which from this point of view are just special values of /3: /3 = |34[ I36| and /3 = 1 [37|. This general 
form of the potential appeared already (in a different parametrization) in Ref . '36' , but only the case 
P — was discussed there. Kimm et al. 35 have studied the general case in flat space and found three 
kinds of flux tube solutions characterized by the asymptotic value of f{r), which (for /(O) — 0) may be 
either 0, /3 or oo. 
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If this Chern-Simons gauged 0(3) model is coupled to gravity, we find that the ill-behaved potential 
from the Higgs system becomes now bounded from below. First we integrate H5.17|l either directly or 
using H5.19|l and (|5.26|) to get: 



^(1*1) = 



27rG3eV 

k(1+/32)2 



A 



(/32 - |<i>|V^2)2 
(l+|$|7Ai2)2 



(5.28) 



where the role of the integration constant is played by A defined by A — (1 (3'^)'^ a / ji^ . Then we find 
from (|5.1U|) a new potential with A as an additional (non-negative) free parameter: 



4 4 
6 7* 



8^*2(1-1-/32)2 



|$|2(/32_|4>|2/^2)2 ^^^^4 
(1-H|$|2/Ai2)4 (I+/32; 



(1 |$|2//i2)2 



(5.29) 



This potential is clearly bounded from below, thus restoring the stability. Moreover, it has always a 
local minimum at |$| = w = /3/i, which breaks the U(l) local symmetry. The extremal points at |<i>| — 
and |<i>| — > 00, which are U(l)-symmetric, may be either minima or maxima depending on the following 
conditions: 



13^{1 + [3^) + 47rG3Ai^(/3* - A) > minimum at 
1 + -I- 47rG'3^^(l - A) > minimum at 



m - 



= 
00 



(5.30) 



The potential of the gravitating Chern-Simons type of the usual Higgs system, is obtained in the limit 
/Lt — > 00 (with A/i** and fifji kept finite). 

Actually a special case of this potential for A = was recently obtained by Abou-Zeid and Samtleben 
021 from a different direction of three dimensional N — 2 supergravity theory. We will see that A 7^ 
yields interesting solutions as well. 

Now we concentrate in rotationally symmetric solutions, i.e. we use the Nielsen-Olesen ansatz 
including the time component of the gauge potential: 



$ = nf{r)e 



171 if 



Afj^dx^" = Ao{r)dt + A{r)d(p 



(5.31) 



with the additional requirement that all the metric components depend on r only. 

The field equations for self-dual rotationally symmetric solutions are quite similar to the ones of 
the GMH model with two additional equations: For Aq and for the metric component L^(r) (Lr now 
vanishes). The dependence of Aq on / is obvious from Ea. (|5.26|l . For L^pir) we have: 



ff2. 



(5.32) 



where the dependence of ^ on / is obtained from H5.28|l . The self-duality equation for the scalar field is 
now [t] = +1 as usual): 

rf 

-J- = eA{r) - n (5.33) 

where for brevity we denote A ^ A — AqL^. The equation for A(r) may be easily obtained from (|5.15|) . 
expressing £ and Aq in terms of /: 



A' 
Wr 



/32_/2 



2«;2(l + /32) 1 + 



27rG3/i^ 



(H-/32)2 



A 



P 



(5.34) 



Einstein equations for the metric field (|5.21|) reduce in this case to the following: 



(logH)'-5^Ao(log/)' 







(5.35) 



The reflection symmetry observed in the GMH system exists in the GCSH system as well, provided A is 
also rescaled according to its /3 dependence. We may therefore limit ourselves to solutions with /(O) — 
as before. Clearly there are three kinds of solutions classified by the three possible values which /(oo) 
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may take, according to the different minima of the potential. However, these three values U{Q), U{v) 
and U (oo) are generally different from each other, and each of them vanishes for a different value of A. 
Since vanishing value of the potential minimum is an essential property of our localized solutions, we 
encounter here a situation which differs from that in flat space |35| . where all three possible boundary 
conditions are realized for the same potential. In the gravitating case, each kind is realized for a different 
value of A. The flux tube solutions with quantized flux and the usual boundary condition /(oo) = /3 
exist for A = only. Representative solutions are depicted in Figs l3l4l The two other solution types 
with /(oo) = 0,00 exist for A — (3^,1, respectively and are depicted in Figs l5l8l We may therefore 
summarize the situation by the following relations: 



A = /3* 
A = 
A = 1 



lim /(r) = 

r — >oo ' 

lim /(r) = P 

r — *oo ' 

lim /(r) = oo 



(5.36) 



A special case is A = /3 = 1 which allows both boundary conditions (/(oo) = 0, oo) for the same 
potential. 

Note that all three kinds of boundary conditions ensure also a vanishing asymptotic value of the 
angular momentum density £, which is a necessary condition for a finite angular momentum. The gauge 
potential Aq need not always vanish asymptotically. It docs only in the /(oo) — (3 case. 

For all these three kinds of boundary conditions, the matter fields tend asymptotically to their 
vacuum values. Thus asymptotically, the geometry of space-time will be rotating conic metric. The 
angular deficit may be easily obtained from (|5.35|l to be: 



6lp 
2^ 



nAo(O) + Ao(oo)(eA(oo) - n)] 



For vanishing c (or ^o(oo) or A, i.e. quantized flux) we find 

5ip SttG^k 



27r 



Ao{0) = 87^G3^I'\n\ / f£i{f)df 



(5.37) 



(5.38) 



that is the same as the result (j4.6|l for the GMH model. For the 0(3) case under consideration here, we 
find the same expression as we had in the case of 0(3)-GMH: 



2^ 



lim f{r) — (3 



(5.39) 



The angular deficit for the other solutions without fiux quantization (i.e. non-vanishing c or A) can be 
easily calculated in the 0(3) model from (|5.37|l : 



6(p 

277 



A^^G3^i^f3^e\A{oo)\ 



1^ = 4.G3M^ 
Ztt 



e^(oo) 



(3^ 



lim f{r) = 

r — *C30 

, lim f{r) = 00 



(5.40) 
(5.41) 



Note the dependence on the unquantized magnetic fiux represented here by ^(00). Thus we have in 
these cases a continuum of values for the angular deficit for a given potential and n and /(oo) values. 



We end this section by computing the angular momentum J of the rotationally symmetric solu- 
tions, which we get by integrating over all two-space. The angular momentum also determines the 
asymptotic value of the metric component L^{r) such that: 



(1 - ^)i^(oo) = -4GJ (5.42) 

ZTT 
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This relation is most easily obtained in a Gaussian normal coordinate system, in which the line element 
is written as: ds^ — —dp^ + hapdx^dx^ {x" = t,(p). Since hap depends only on p, the components of 
the Ricci tensor have the following simple form in terms of the extrinsic curvature K^: 



It turns out that we can compute the angular momentum directly without even the self-duality assump- 
tion, due to the simple identity which holds for solutions of the flux tube form (|5.31|l : 

eT^ = in- eAir))f (5.44) 

Now we use the time component of the Chern-Simons equation (j5.3|l . to write in terms of B and get: 

J=— (n-eA(r))A'dr^—(2n-eA(oo))A(oo) (5.45) 

Solutions which break the U(l) symmetry asymptotically have quantized angular momentum of: 



TTK 



J = -^n^ (5.46) 



e 



while otherwise, there is no flux quantization and we cannot say about the angular momentum more 
than the right hand side of H5.45|l . Note that the angular momentum in (|5.4t)|l is a very general result 
independent of all details of the theory, except the existence of symmetry breaking vacuum in the 
potential. We note further that all solutions have the same sign of angular momentum, which is a 
manifestation of parity violation. 



6 Correspondence between Maxwell and Chern-Simons Type 
of the Generalized Higgs Model 

We note that in flat space, the electric energy density of the Chern-Simons field plays an equivalent role 
to that of the magnetic energy density in the GMH model. Moreover, due to H5.15|l . there is a relation 
between the magnetic and electric terms, thus enabling us to eliminate the scalar (electric) potential 
from the equations and to give them a form identical to those of the GMH model. 

Comparison of the Bogomolny constraints in both cases, equations H2.17|l and (|5.18|) . gives a relation 
between the magnetic field in the GMH model and the time component of the gauge potential in the 
GCSH model: 

£2B = nAo (6.1) 

Further comparison of the magnetic field appearing in this equation with the one in (|5.15|) (with ^ = 0) 
or the potentials f H2.11|l and H5.10|l (with £ = 0)), yields the following characteristic relation for the 
generalized Maxwell-Higgs model: 



K 



2 



Thus, all the flat space self-dual solutions to the D — 3 system, studied in section 5, are also self-dual 
solutions to the D = 4 GMH model of section^ provided we use £2 consistent with H6.2|) . Consequently, 
the curves in figures 1^15171 represent also the corresponding fields of the GMH model with the triple well 
potential (|3.13() and the appropriate dielectric function mentioned together with that potential. 

However, the presence of an electric field in the D = 3 Chern-Simons theory, in addition to the 
magnetic field, is responsible for the existence of angular momentum, which calls for corrections to this 
simple correspondence when gravity is considered dynamical. It is straightforward (although tedious) to 
show that there does not exist spinning gravitating self-dual flux tubes in the GMH model. Therefore, 
we have only static versus stationary correspondence, and all we can hope for is equivalence of the 
two-dimensional spatial metrics in addition to the scalar and vector correspondence. This turns out to 
be the case, provided we also impose c = a = X = 0. 
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Equation (|6.1(l is evidently still valid in the gravitating case and should be accompanied hy B — B 
where B is the GMH magnetic field defined by H2.9|l and B - in ()5.8|l . Thus we get the relation instead 
of (lOl: 



^^=e^|<i>|^£i(|$|) + Kf (6.3) 
Note that now the potentials are not equal in both cases, but there is an ^-dependent difference: 

Um = Ucs - (6-4) 

As for the metric tensor, it is easy to see that both expressions for the Ricci scalar, equations H2.27|l 
and H5.2U|I . become identical due to equation Ht).l|l . thus resulting in equal 2-metrics in the two cases. 

Actually, those ^-dependent modifications should be expressed in terms of functions that appear in 
the Lagrangian, and should be considered as functions of |$|. We will not do this explicitly, since it is 
straightforward and does not add any more insight into the physical picture. 



7 Topological Charge 

The quantization of magnetic flux is intimately connected with the fact that the minimum of the 
potential breaks the local U(l) symmetry. Consequently, the magnetic flux number is the index, or 
winding number, of the map defined by the scalar fields from infinite distance from the fiux tube to the 
vacuum manifold. 

If however, the target space is compact, there exists a further possibility of homotopy classification 
of the maps defined by the fields from "real" space to target space. We can use, as a ground of the 
present discussion, the well-known notions from sigma models defined in three spacetime dimensions 
|14j . The GCSH model is naturally three dimensional, but the GMH model should be thought of here 
as three dimensional as well. In this case, the scalar fields map the two dimensional physical space into 
target space. We have therefore an identically conserved topological current defined by: 

where Or is the target space volume, and .^d^l) is defined by: 

This is a gauge invariant generalization of the standard non-gauged sigma model current, whose time 
component integral over all space is the winding number, or index, of the map defined by the sigma 
model fields. Note that T{\^)\s connected to ^(|$|) (defined in Ea. (|T^ ) by eJ^(|$|) = 
up to a possible additive constant. By the field equations of the GCSH model, it is also related to Aq 
by kAo(|$|) = — 77eJ^(|$|), up to an additive constant. Unlike previous authors 30, 35, 36jj37j, we fix 
the additive constant of J- such that J- vanishes as j^j — ^ oo, irrespective of the boundary conditions 
imposed on the solutions. This way, we can treat uniformly all sigma models differing only by the 
potential term. For the special 0(3) "dielectric function", given in (|2.22l) . we may write the topological 
charge (and current) in terms of ^id^j) only, since JF(|$|) — — /i^-\/£i(|<I>|)/2. In this case, we have also 
f^T = TT/i^ and the apparent dependence on the scale /i disappears. 

The fact that this gauge invariant current has still a topological meaning is clearly seen by the fact 
that it is actually the non-gauged sigma model current with an addition, which is a divergence of an 
anti-symmetric tensor: 

K>^ ^ --f^ {^£lm){^,'^>n^.^) + 2e\/,{Tm)A,)) (7.3) 
2i2TVl.9l 

Indeed, the first term of (|7.3|1 can also be written as a divergence of an anti-symmetric tensor, so we 
have the alternative form: 

2nTxAQ\ \\^\ 



K" = V^Kf"", Kf"" = ( j^i^*dx'S> - $5a$*) + "^eAx ] (7.4) 
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Integration of the time component of the topological current provides us with the topological charge 
T, which characterizes the solutions from the homotopy point of view. For our stationary solutions, we 
write the two equivalent expressions: 

T = j d?x (ze'^^fi(|$|)(A$)* - 2e^\Bnm 



J d'x^=E,m){dr^r{d,'^) - ^H^) f^^^ Adx^ (7.5) 

where w is the asymptotic value of |$|. We stress the difference between w, which corresponds to a 
property of a solution, and v which is a parameter in the potential, w need not be equal to v (although 
it may be), as happens for solutions of the GCSH model discussed in section|Sl or solutions to the GMH 
model with the potential H3.13|l . which may have all three asymptotic values w = 0,v, oo. 

Note that T cannot be expressed as a surface integral (actually line integral in 13 — 2 + 1) since 
the winding number, i.e. the (non-gauged) sigma model contribution of K^'^, does not satisfy Stokes 
theorem. The second term is however a boundary term, and we note that it vanishes for solutions for 
which |<I>| ^ oo as r ^ oo. In the 0(3) model, these are solutions which do not break the U(l) symmetry 
asymptotically. Another possibility is solutions which do break the U(l) symmetry asymptotically, and 
thus do not cover all target space. This kind of solutions have non-integer winding number, but the 
topological charge has also a magnetic flux contribution, which may compensate for this fact. 

This is exactly what happens for the self-dual cylindrically symmetric solutions of the 0(3) GMH 
model, where we find that in all cases T = n. If the solution tends to a symmetry breaking vacuum, 
the boundary term of H7.5|l does not vanish and evaluates to n/(l + uj'^), where we define uj — w/ fi. The 
winding number, which is the first term, is easily found to be nw^/ (1 -I- ux^). They both add up to give 
T = n. In the case ui — > oo, where the solution does not break asymptotically the U(l) symmetry, the 
flux contribution vanishes but the first is an integer and we get T — n again. The third possibility, a; = 
is however an exception: In this case, the winding number vanishes but we have a flux contribution to 
get T = eA{oo). 

Another property, which we can note here, is the relation with the Euler number of the 2-surface 
generated by the self-dual GMH solutions. The Euler number is given by: 

X - / <fxVW\Rh) (7.6) 

while for (say) self-dual solutions (77 = 1), we may write the topological charge as: 

T = '^ I d'x^\{£,imY'iD.<i>nD,<i>) + 2£2im)B') - -^Hv) <f^ ^Adx^ (7.7) 

By the field equations H2.11|l and H2.14|l . we see that the first term is just proportional to the Euler 
number, so we have a simple relation between the topological charge, the Euler number and the magnetic 
flux, which we denote here by 5": 

For the 0(3) system, (|7.8|l yields the following relation: 

2'^^+2i^+Ti^^^o ^'-'^ 

For asymptotically conic space, the Euler number is related to the deficit angle x — 5lp/2t:^ so 
Ea. (|7.8|l is easily verified by using H4.7|l . 

Next we turn to the GCSH model, where things are quite similar. The topological current and charge 
are still given by Eas. (|7.1|) - (|7.5|) . However, since the self-duality condition is now modified by replacing 
Ai by Ai, it is useful to write the topological charge also in a way that is ready for direct use of the 
modified self-duality. We therefore write the topological charge density as: 

= - (z£i(|a>|)(A$)*(A$) + 2ed,{Tm)AoL,)) + ^' (7.10) 

2ilTy'\g\ "T 
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Note however that the surface term vanishes only for quantized flux sohitions (c — X — 0, Aq{w) = 0). 
If we proceed along similar lines as for the GMH case, we easily find the same relation between the 
topological charge, Euler number and the magnetic flux, i.e. Ea. H7.8() . It simplifies again for the 0(3) 
model to (171^1 . 

Finally, we use this relation in order to calculate the topological charge for the solutions with quan- 
tized as well as non-quantized flux, which is obtained directly from Ea. (|7.5(l : 

{eA{oo) , w = 
n ^ w = v (7.11) 

n , w ^ oo 

A nice check of our results is that we can easily reproduce by substitution in 1)7.9(1 the expressions 
l)5.39|l - H5.41|l for the deficit angles. 



8 Outlook 

We have shown that the inclusion of dielectric functions gives the possibility of dramatic changes in 
the long range behavior of scalar, gauge and metric fields around a cosmic string or vortex. This was 
obtained even for very simple and symmetric models, like 0(3) with the potential (|3.13ll . We believe 
that the astrophysical and cosmological consequences deserve further studies. 

Another continuation and generalization of our work is to consider theories with a more extended 
field content. One natural possibility is to gauge a U(l)xU(l) subgroup of the global symmetry group of 
a CP(2) non-linear sigma model noticing that CP(2) can be parametrized by two complex coordinates. 
One can therefore expect that in analogy with the way the 0(3) model (with ~ CP(1) as target 
space) generalizes the ordinary cosmic string, the gauged CP (2) model generalizes the superconducting 
cosmic string |43l 144) . 
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Figure 1: The solution to Eas. H3.3|I . H3.4|l . for £2 = 1- B is the angular coordinate on S^. We also 
show the dimensionless magnetic field. The metric field, which is constant in this case, is included for 
comparison with the following figures. The parameters used are n = — 1, [3=1. The dimensionless 
length coordinate is defined by x = efir. 
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Figure 2: The solution to Eas.H4.1 |) - (|4.3|) . for £2 = 1- The parameters used are n — —1, (3=1, SirGfi'^ 
0.5. Compare with Figure 1. 
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Figure 3: The solution to Eas. H5.32|) - (|5.35|l . with the boundary condition /(oo) = (3. The parameters 
used are n = —1, /3 = 1, n/e^— 1, SnG^fi'^ = 0, i.e. gravity is not included. Notice that we magnify the 
magnetic field by a factor 10 here and in the subsequent figures. The gauge field tends to its asymptotic 
value slower than in the previous case as is seen from the value of eA at x = 20 which is somewhat 
above -1. The curves in this figure represent also the corresponding solution for the GMH system with 
the triple well potential - see remark below Eq. (|(j.2p . 
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Figure 5: The solution to Eas. H5.32l) - (|5.35l) . with the boundary condition /(oo) = oo. The parameters 
used are n = —1, /? = 1, K/e/i — 1, SirGsfj,^ — 0, i.e. gravity is not included. The curves in this figure 
represent also the corresponding solution for the GMH system with the triple well potential - see remark 
below Eq. |jn21l- 
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Figure 6: The solution to Eas. H5.32(l - (|5.35(l . with the boundary condition f{oo) = oo. The parameters 
used are n = —1, (3—1, K/e/i = 1, SnGsfi^ = 0.0625, i.e. gravity is included. Compare with Figure 5. 
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